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ABSTRACT 


Both  the  total  electron  content  of  the  ionosphere  (TEC)  and 
the  critical  frequency  of  the  F2  layer  (foF2)  exhibit  large  day- 
to-day  variations  during  quiet  and  active  geomagnetic  periods.  It 
is  of  great  interest  to  ascertain  whether  good  correlation  exists 
between  TEC  daily  variability  about  the  monthly  mean  and  foF2 
variations.  With  the  availability  of  the  global  GPS  constellation 
to  provide  instantaneous  time-delay  values  such  a  correlation  may 
enable  the  improvement  of  HF  short-term  predictions  using  passive 
monitoring  of  TEC. 

To  determine  the  correlation  one  year  of  TEC  and  foF2  data 
were  collected  from  June  1995  until  May  1996.  The  TEC  data  was 
determined  from  GPS  time-delay  measurements  at  Matera,  Italy, 
obtained  from  the  data  base  of  the  International  GPS  Services  for 
Geodynamics.  The  foF2  measurements  came  from  Rome,  Italy. 

The  analysis  showed,  that  for  large  percentages  of  the  time 
very  good  correlation  exists  between  TEC  and  foF2  short-term 
variations.  The  correlation  coefficient  varies  from  0.7  -  0.8 
during  the  summer  months  to  about  0.5  -  0.6  during  the  winter.  A 
study  of  the  diurnal  dependence  of  the  correlation  indicates  that 
better  correlation  exists  during  day-time  than  night-time. 

The  high  correlation  between  TEC  and  foF2  indicates  that 
real-time  ionospheric  HF  prediction  improvements  are  feasible 
when  using  transionospheric  time-delay  measurements. 
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1.  INTRODUCTION 

HF  radio  communication  depends  on  the  ability  of  . the 
ionosphere  to  return  the  radio  signal  incident  on  it  back  to 
Earth.  Predictions  of  ionization  levels  in  the  various  ionosphe¬ 
ric  regions  are  derived  from  models  and  are  used  as  a  basis  for 
planning  and  frequency  management  of  HF  radio  systems  worldwide. 
Uncertainties  or  inaccuracies  in  the  models  of  the  ionosphere 
have  long  been  known  to  be  one  of  the  major  causes  if  not  the 
major  cause,  for  inaccuracies  in  the  calculated  propagation 
characteristics.  This  is  particularly  true  for  those  applications 
of  ionospheric  predictions  involving  timescales  that  are  less 
than  the  monthly  average  or  monthly  median. 

To  reduce  average  monthly  RMS  errors  in  predictions,  adaptive 
techniques  that  use  real-time  observations  to  correct  model 
biases  have  been  devised.  Examples  of  such  are  (1)  sounding  the 
desired  communication  path  prior  to  information  transmittal  or 
(2)  sounding  the  vertical  ionosphere  to  determine  a  real-time 
model  reference  point  and  then  to  adjust  the  model  for  the 
future,  until  the  next  reference  sounding.  While  the  first  is 
appropriate  for  fixed  point  communication  paths,  it  is  cumbersome 
and  causes  EM  interference  over  wide  geographic  areas,  whereas 
the  second  technique,  while  not  causing  interference,  may  not 
yield  values  representative  of  ionospheric  conditions  at  the 
reflection  point  of  an  oblique  path. 

A  method  that  is  potentially  global  in  nature  involves  the 
monitoring  of  satellite-emitted  signals  which  yield  information 
on  the  ionospheric  parameters  along  the  propagation  path  such  as 
time-delay  or  TEC  and  converting  such  information  into  the  HF 
propagation  parameters  of  interest  (foF2).  The  advantage  of 
monitoring  satellite-emitted  signals  is  the  fact  that  it  is 
passive  for  the  potential  user,  and  the  existence  of  a  global 
network  of  satellites  (e.g.,  the  Global  Positioning  System,  or 
GPS)  affords  the  possibility  of  global  coverage. 
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To  assess  the  possible  immprovement  of  HF  short-term 
predictions  from  transionospheric  measurements,  a  correlation 
study  between  TEC  daily  variability  about  the  monthly  mean  and 
f oF2  variability  was  conducted.  Several  months  of  TEC  data  taken 
in  Haifa,  Israel  during  1980,  as  well  as  GPS  time-delay 
measurements  taken  during  the  summer  of  1992  in  Jerusalem,  Israel 
were  correlated  with  foF2  measurements  from  Cape  Zevgari,  Cyprus. 
The  results  of  this  study  were  reported  and  published  in  various 
scientific  conferences  [1],[2],[3],[4],[5],  scientific  journals 

[6]  and  the  Final  Technical  Report  of  contract  DAJA-93-C-0035 

[7] . 


The  above  results  showed  that  very  good  correlation  exists 
between  TEC  and  foF2,  indicating  that  real-time  HF  prediction 
improvements  by  passive  means  are  feasible.  However  the  previous 
studies  were  based  on  a  limited  amount  of  data  and  thus  it  was 
decided  to  increase  the  data  base  by  collecting  new  foF2  and  TEC 
measurements  for  at  least  one  more  year.  In  this  report  we 
present  the  results  of  the  new  study,  which  was  conducted  during 
the  minimum  phase  of  the  last  sunspot  cycle. 

2.  EXPERIMENTAL  DATA 


The  experimental  data  is  divided  into  two  periods:  1994-1995 
and  1995-1996.  The  first  were  collected  from  June  1994  until  May 
1995.  The  TEC  data  were  obtained  from  GPS  observations  taken  at 
the  National  Physical  Laboratory  in  Jerusalem,  while  the  foF2 
measurements  were  from  Nicosia,  Cyprus  and  were  received  from 
World  Data  Center  A  in  Boulder,  Colorado.  Unfortunately  most  of 
the  data  received  both  from  Jerusalem  (for  the  period  July- 
November  1994)  as  well  as  from  Boulder  (from  September  1994  until 
May  1995)  were  either  unreliable,  missing  or  incomplete  and  could 
not  be  used  for  the  correlation  analysis. 

When  it  was  realized  that  we  would  not  be  able  to  use  the 
Cyprus  data,  it  was  decided  to  use  other  sources.  GPS  time-delay 
measurements  can  be  obtained  from  the  data  base  of  the 
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international  GPS  services  for  Geodynamics  (IGS).  The  IGS  data 
base  contains  measurements  from  some  50  locations  around  the 
world  and  can  be  accessed  via  the  internet.  The  nearest  location 
to  Israel  is  in  Italy  and  thus  from  June  1995  until  June  1996  we 
obtained  GPS  time-delay  data  from  Matera,  Italy  (40.65  N  16.7  E) 
while  the  foF2  measurements  came  from  Rome,  Italy  (41.9  N  12.5  E) 

The  TEC  were  determined  from  the  raw  IGS  data  with  a  software 
package  developed  by  the  Philips  Laboratory,  Hanscom  Field,  Mass 
and  adapted  for  use  on  the  Technion  computer  [8]. 

3.  CORRELATION  RESULTS 

For  each  month  of  observations,  the  hourly  values  of  the 
variability  of  both  the  TEC  and  foF2  were  determined.  The 
variability  is  calculated  by  substracting  the  monthly  average 
value  from  each  hourly  value  and  dividing  by  the  monthly  average 
value.  Another  parameter  of  interest  is  the  slab  thickness,  which 
is  proportional  to  the  ratio  TEC/(foF2)2.  The  slab  thickness  can 
be  calculated  for  each  hour  of  observations  and  its  variability 
determined  in  the  same  way  as  for  TEC  or  foF2. 

Both  the  variation  in  foF2  and  TEC  are  noisy  on  an  hour-to- 
hour  basis,  which  is  caused  by  measurements  and  data  reduction 
errors  rather  than  by  physical  phenomena.  To  reduce  the  noise, 
the  variations  in  TEC  and  foF2  were  smoothed  by  a  3  hour  running 
mean.  Crosscorrelation  analysis  was  performed  on  the  foF2  and  TEC 
variabilities  both  for  the  unsmoothed  and  smoothed  data. 
Autocorrelation  analysis  was  also  performed  on  the  unsmoothed  TEC 
and  foF2  data. 

A  complete  set  of  results  for  the  period  June  1995-May  1996 
are  given  in  Appendix  A.  The  following  figures  are  shown  for 
every  month: 

a.  Hourly  values  of  the  variability  of  foF2  and  TEC. 

b.  Smoothed  hourly  values  of  the  variability  of  foF2  and  TEC 
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c.  Cross-correlation  function  between  hourly  values  of  foF2 
and  TEC  for  the  unsmoothed  data. 

d.  Auto-correlation  functions  for  the  unsmoothed  foF2  and 
TEC. 

e.  Hourly  values  of  the  variability  of  the  slab-thickness  r. 

The  cross-correlation  coefficients  for  both  the  smoothed  and 
unsmoothed  data  for  the  period  June  1995-May  1996  are  given  in 
Table  1.  It  can  be  seen  that  the  correlation  coefficients  for  the 
unsmoothed  data  vary  between  0.78  for  June  1995  to  0.53  for 

February  1996.  The  correlation  coefficient  for  the  smoothed  data 
vary  between  0.84  and  0.58  for  the  same  months. 

The  diurnal  variations  of  the  correlation  coefficients  for 

the  various  seasons  are  shown  in  Figure  1 .  Each  point  in  the 
figure  represents  correlation  results  for  a  4  hour  block.  It  can 
be  seen  that  the  diurnal  variations  are  similar  for  all  the 
seasons  besides  the  winter  of  1996.  Minimum  correlation 

coefficients  accur  at  night  near  sun  rise  and  a  broad  maximum  is 
observed  during  the  day-time.  There  is  no  apparent  diurnal 

variation  for  the  winter  months.  The  reason  for  this  is  unclear. 
The  highest  correlation  coefficients  for  both  day-time  and  night¬ 
time  are  observed  in  the  summer  of  1995. 

4.  CONCLUSIONS 

The  high  cross-correlation  coefficient  for  foF2  and  TEC 
raises  the  possibility  that  real-time  TEC  measurements  may  be 
used  to  update  foF2  value  determinations.  The  cross-correlation 
may  even  be  higher  if  the  geographic  subionospheric  point  of  the 
TEC  measurements  is  closer  to  the  geographic  point  of  the  foF2 
measurement,  which  introduces  an  error,  in  addition  to  the 
possible  inherent  measurements  un-certainties.  TEC  measurements 
utilizing  satellite  emitted  signals  are  passive  in  nature  and  do 
not  burden  the  electromagnetic  spectrum.  In  addition,  the  availa¬ 
bility  of  the  global  GPS  constellation  to  provide  instantaneous 
time-delay,  or  eguivalently  TEC,  values  could  provide  an  instan- 
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Table  1 


TTALY  1995/6 


MONTH 

CORR  COEF 

NO  SMOOTHING 

CORR  COEF 

3 -POINT  SMOOTHING 

JUN  1995 

0.78 

0.84 

JUL 

0.75 

0.82 

AUG 

0.70 

0.76 

SEP 

0.74 

0.80 

OCT 

0.59 

0.67 

NOV 

0.64 

0.75 

DEC 

0.66 

0.76 

JAN  1996 

0.63 

0.74 

FEB 

0.53 

0.58 

MAR 

0.66 

0.72 

APR 

0.63 

0.68 

MAY 

0.60 

0.64 

CORR  COEF 


Figure  1 
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taneous  updating  of  foF2  models  on  a  global  basis  as  well  as  on  a 
regional  basis.  Such  capability  is  important  for  HF  communication 
along  short ,  medium  and  long  range  paths . 
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APPENDIX  A 

Experimental  results  for  the  period  June  1995-May  1996. 

The  following  figures  are  shown  for  every  month: 

a.  Hourly  values  of  the  variability  of  foF2  and  TEC 

b.  Smoothed  hourly  values  of  the  variability  of  foF2  and 
TEC. 

c.  Cross-correlation  function  between  hourly  values  of 
f oF2  and  TEC  for  the  unsmoothed  data. 

d.  Auto-correlation  functions  for  the  unsmoothed  foF2  and 
TEC  variations. 

e.  Hourly  values  of  the  variability  of  the  slabthickness . 
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1-31  August  1995 
3  point  smoothing  R  =  0.76 
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1-30  September  1995 
3  point  smoothing  R  =  0.80 
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1-31  October  1995 
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